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The Bose Hubbard model (BHM) is an archetypal quan- 
tum lattice system exhibiting a quantum phase transi- 
tion between its superfluid (SF) and Mott-insulator (MI) 
phase. Unlike in higher dimensions the phase diagram of 
the BHM in one dimension possesses regions in which 
increasing the hopping amplitude can result in a transi- 
tion from MI to SF and then back to a MI. This type of 
re-entrance is well known in classical systems like liq- 
uid crystals yet its origin in quantum systems is still not 
well understood. Moreover, this unusual re-entrant char- 
acter of the BHM is not easily captured in approximate 
analytical or numerical calculations. 



1 Introduction Re-entrance is a novel phenomenon 
in which a succession of transitions between two phases 
A and B, such as A-B-A-B, can occur by monotonically 
increasing just one parameter. Such a series of transitions 
implies a non-linearity which allows a single parameter 
to produce opposite effects depending on its initial value 
and is often a reflection of the interacting character of the 
system. In the context of classical thermal phase transi- 
tions the natural parameter to vary is temperature. It is 
then expected that the low temperature phase A will be or- 
dered, while increasing the temperature will drive the sys- 
tem to a disordered phase B. However, the appearance of 
re-entrance means that increasing the temperature can in 
fact unexpectedly stabilize the ordered phase A again. This 
sequence of phase transitions has been observed in liquid 
crystals between the A = Smectic (ordered) phase and B = 



Here we study in detail the predictions of three differ- 
ent and widely used approximations; a multi-site mean- 
field decoupling, a finite-sized cluster calculation, and 
a real-space renormalization group (RG) approach. It is 
found that mean-field calculations do not reproduce re- 
entrance while finite-sized clusters display a precursor to 
re-entrance. Here we show for the first time that RG does 
capture the re-entrant feature and constitutes one of the 
simplest approximation able to do so. The differing abil- 
ities of these approaches reveals the importance of de- 
scribing short-ranged correlations relevant to the kinetic 
energy of a MI in a particle-number symmetric way. 
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Nematic (disorded) phase with increasing temperature [Q] 
[2], and similar behaviour also occurs in frustrated classical 
spin systems J3). 

The appearance of re-entrance in the zero-temperature 
phase diagram of a quantum lattice system is much less 
well understood. In this work, we study the Bose-Hubbard 
model (BHM) which is a minimal many-body Hamiltonian 
where the competition between the kinetic and repulsive 
interaction effects gives rise to a quantum phase transi- 
tion [4)- For weak interactions the bosons are completely 
delocalized leading to a superfluid (SF) phase, while for 
sufficiently strong interactions, and a commensurate den- 
sity, the bosons become localized and enter the Mott insu- 
lator (MI) phase. The BHM has long been used to describe 
Josephson junction arrays in which several superconduct- 
ing island are coupled together Q. More recently it has 
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received considerable attention and experimental promi- 
nence due to its near perfect realization with cold atoms 
trapped in optical lattices [6,7,8,9 10 L 11]. 

The qualitative structure of the BHM phase diagram, 
such as the existence of Ml lobes depicted in Fig. Q2 a X 
was worked out some time ago by Fisher et al. Ifl2l . They 
showed that the tip of the lobes undergo a transition be- 
longing to the universality class of the XY model in d + 1 
dimensions, implying that in one dimensions (ID) it is a 
Kosterlitz-Thouless (KT) lfl3ll transition. Only much later 
was it discovered that the Ml lobes in ID actually dis- 
plays re-entrance at low filling factors. Specifically at some 
constant chemical potential, near the tip of the lobe, a re- 
entrant sequence of zero-temperature quantum phase tran- 
sitions occur, with A = MI and B = SF, as the coherent 
hopping amplitude is increased from zero. This is again 
surprising since it demonstrates that increasing the hopping 
amplitude, which in isolation favours the itinerancy of the 
bosons, can instead favour their localization under certain 
circumstances. This delicate behaviour is not easily cap- 
tured in commonly used approximations, a fact connected 
to the known difficulty of describing the tips KT transi- 
tion. Instead the existence of re-entrance was firmly es- 
tablished using sophisticated numerical tools such as den- 
sity matrix renormalization group (DMRG) II14II151 . high- 
order strong-coupling perturbation (SCP) expansions fTSl 
[TTll . and quantum Monte Carlo (QMC) simulations 11811191 
l20l . However, an understanding of the underlying mecha- 
nism leading to re-entrance has yet to be achieved. Our aim 
here is to glean insight into the key physical properties re- 
quired for its emergence in the BHM in ID by analyzing 
the capabilities of relatively simple and intuitive approxi- 
mation methods to capture it. 

The structure of this paper is as follows. In Sec. [2] we 
discuss the essential properties of the MI and SF phases 
and employ several qualitative arguments to understand 
the structure of the BHM phase diagram. This is followed 
by a discussion of the exact form of phase diagram for 
the BHM in ID, displaying re-entrance, as determined by 
earlier DMRG, SCP expansions, and QMC calculations. 
Then in Sec. [3] three different widely used approximation 
methods are applied to the BHM. This includes a multi- 
site mean-field decoupling l2T]|22ll23ll4ll24ll25l. a finite- 
sized cluster calculation [26 , 27 , 28 , 29 1 and a renormaliza- 
tion group (RG) approach BTI . Finally, in Sec. [4] several 
conclusions are summarized from the capabilities of these 
methods at describing re-entrance. 



2 Bose-Hubbard model The BHM describes a 
bosonic system in a regular lattice where particles can 
tunnel to nearest neighbors and interactions occurs be- 
tween bosons at the same lattice site. The Hamiltonian in 



the grand canonical ensemble is (taking fi = 1) 

# = -*£ + h - c -) - 1* J2 h i 

+^X>j(^-i), (i) 

where bj , are the bosonic creation and annihilation oper- 
ators at site j and fij = b^bj is the corresponding number 
operator. Here t > is the hopping amplitude of bosons 
between neighbouring sites, /i is the chemical potential and 
U > is the repulsive on-site interaction strength. In the 
non-interacting limit, U = 0, the Hamiltonian Eq. Q] is 
diagonal in the basis of single-particle plane waves with 
a tight-binding dispersion relation e(fc) = — 2icos(fca), 
where a is the lattice spacing. In this case, the ground- 
state is characterized by all the particles condensing into 
the k — state. The number of particles in the ground- 
state diverges whenever [i < —2t. In the following we fix 
our energy scale by U and consider the phase diagram in 
the (t, (j,) plane where for simplicity we label the dimen- 
sionless ratios as t/U —> t and fi/U —tfi,. 

In the t = limit, the energy eigenstates possess a well 
defined occupation at each lattice site. The ground state is 
uniformly filled with n bosons in each site giving a com- 
mensurate density so long as the chemical potential ful- 
fills El 

n — 1 < < n. (2) 

In the case fi > 0, the ground state at zero hopping is filled 
with at least one boson per site. Excitations above this state 
can be constructed by adding (particle excitations) or re- 
moving (hole excitations) a boson to the ground state at any 
site. The energy of these lowest particle and hole excitation 
with respect to the ground state energy is A p (fi) = — fi+n 
and Ah (fx) = \i — n + 1. For a non-integer fi, we see from 
Eq. © that A p (fj,) > and A h (fj,) > 0. Thus, a defin- 
ing characteristic of the MI phase is the gap induced by the 
interaction for both particles and hole excitations ll32l . 

Keeping fi > constant, the particle and hole gap de- 
creases when hopping is increased from the t = limit. 
Depending on /j, there is a hopping value at which one 
of the particle or hole gaps becomes zero. At this point, 
the system undergoes a phase transition to the gapless SF 
phase. Commonly the system goes from a commensurate 
to an incommensurate density when generically crossing 
the MI-SF phase transition. We refer to the phase bound- 
ary in which the particle gap closes as the particle bound- 
ary and as the hole boundary when the hole gap closes. 
Crossing the phase boundaries from the MI to the SF re- 
gion produces an increase (particle boundary) or decrease 
(hole boundary) in the density. 

In ID quantum fluctuations destroy long-range order 
in the SF phase. Instead, this phase exhibits quasi-long 
range order with algebraic decaying correlations, such as 
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Figure 1 (a) The mean-field phase (t, p) diagram of the BHM in one dimension [21 22 23]. The shaded regions are the 
MI lobes with integer fillings n. These parabolic shaped lobes are surrounded by the SF phase. Some density contours in 
the SF phase are shown to illustrate their tendency to have a negative slope. Note also that the MI lobes are particle-hole 
asymmetric. The hole boundary lower side of the lobe while the particle boundary is upper side, (b) Focusing on the n = 1 
MI lobe the DMRG results from Ref. |[T5l are shown as o, along with two different sets of QMC data from Ref. |fl9l as x 
and Ref. lEUl as O. The solid lines delineate the region found to have unit density using a 12th order SCP expansions ifTTl 
and the large solid point marks the KT tip found from DMRG lfT31 . 



(%b x ) ~ x~ K l 2 with K = Tr/y/p s Kt, where p s is the su- 
perfluid density ll33ll and k = dp/ dp is the compressibil- 
ity l33l . This lack of genuine long-range order shows that 
the SF phase in ID is not a Bose-Einstein condensate. Nev- 
ertheless, it is the non-vanishing p s that defines this phase 
as a genuine SF (34). Within the MI phase correlations de- 
cay exponentially fast with a correlation length £, such as 
(bob x ) ~ cxp(— x/£). These correlations indicate that the 
bosons in a MI state are localized over a typical region of 
length £ contrasting to the SF regime where particles can 
move over the whole system. 

In the SF region, the constant-density region are curves 
in the (p, t) plane due to the need for positive compressibil- 
ity, k > 0. Thus, the particle and hole boundaries surround- 
ing each of the MI regions must collapse in a constant- 
density curve in the SF phase resulting in the well known 
MI lobe structure RTTPHl . It also implies that the MI to SF 
transition at the tip of each of the MI lobes then occurs at 
constant density coinciding with the merging of the parti- 
cle and hole boundaries. The scaling theory developed by 
Fisher et al. Ifl2l showed that generic transitions are in the 
universality class of the Gaussian model, while the transi- 
tion occurring at the tips belongs to the universality class 
of the X F-model in d + 1 dimension. Thus, a KT critical 
point appears at the tips of each MI lobe in ID problem 
so both particle and hole gap close exponentially slowly 
as A(t) = Aexpi-B/y/tKT - t), with A,B real con- 
stants Ifl3l . 

A qualitative description of the phase diagram for the 
ID BHM can be obtained using a single-site MF decou- 
pling approximation [21 , 22, 23 4|. This results in the phase 



diagram shown in Fig. [TJa). The light-blue area indicates 
the commensurately filled MI lobes while the white denote 
SF regions surrounding them. Also shown are some den- 
sity contours in the SF phase. Since in the non-interacting 
limit, equivalent to t —> oo, the particle number diverges 
for any p, this implies that the density contours in the SF 
must have a negative slope whenever the hopping is large 
enough. As can be seen in Fig. Q2a) this feature is repro- 
duced by the MF approximation. This kinetic effect occur- 
ring at moderate values of hopping can be viewed as a pre- 
cursor to re-entrance lfT5l . as we shall see shortly. However, 
this simple single-site MF approximation, which is essen- 
tially equivalent to the infinite dimensional limit, fails to 
predict re-entrant MI lobes. 

3 Capturing re-entrance Re-entrance was not pre- 
dicted for the ID BHM until more sophisticated and accu- 
rate approximations were applied. In particular it was first 
revealed by using 12th order SCP expansions 1171 . After 
this discovery further calculations with DMRG [15| and 
QMC methods 01911201 convincingly established the phe- 
nomenon. Their original results are replotted in Fig. |TJb) 
showing the most accurate, and essentially exact, result for 
the ID BHM phase diagram. A clear consensus from these 
calculations is that the MI lobe in ID and for small filling 
factor n ~ 1, has a triangular pointed shape which is quite 
different from the rounded shape revealed in the MF ap- 
proximation. Moreover the tip of the lobe, also highlighted 
in Fig. |TJb) from the DMRG calculation ITT51 . occurs at a 
much larger value of t than MF predicts owing to the KT 
nature of its transition in ID. Roughly speaking re-entrance 
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occurs due to a combination of this elongation of the tip 
with the kinetic energy driven tendency for all density con- 
tours to eventually have a negative slope. With this view it 
is therefore also natural that re-entrance only occurs in ID 
because in higher dimension the transition at the tip of the 
MI lobe is not a KT critical point and instead occurs at val- 
ues of t increasingly close to the MF prediction. Ultimately 
these two effects cause the hole boundary to become con- 
cave before it intersects with the particle boundary. This 
then makes it possible to traverse, by increasing t alone, 
a series of phases transitions MI-SF-MI-SF at constant fi 
in the region just above the tip. That re-entrance occurs in 
this region, where the interplay between hopping and inter- 
actions are most significant, illustrates how strong correla- 
tions are key to its appearance. We note also that increasing 
the filling n causes bosonic enhancement of the kinetic en- 
ergy and consequently the MI lobe shrinks as « 1 /n. Since 
the tip then occurs at much smaller values of t re-entrance 
disappears for n > 1. For this reason we focus entirely on 
n = 1 filling in this work. 

The sophisticated methods employed to determine this 
structure however, obscure the underlying physical intu- 
ition behind re-entrance owing to their complexity. Our 
main aim here is to capture re-entrance within the simplest 
possible approximation in order to gain further insight into 
the origin of this remarkable phenomenon. We now exam- 
ine in turn some common approaches. 

3.1 Multi-site mean-field decoupling The MF de- 
coupling approach is often a highly successful means of 
determining the qualitative structure of the phase diagram 
for quantum lattice systems. For the BHM the most effec- 
tive strategy is to decouple the hopping term via the factor- 
ization EOS) 

+ (3) 

where (■) denotes the ground state average. Applying this 
factorization leaves the interaction term untouched but re- 
duces the BHM to a sum of identical single site Hamilto- 
nians 

ifi = -2t(0P + <P*b) + 2t\<P\ 2 + ^h(n - 1) - fih. 

It is known that the exact solution of the SF in ID does not 
posses genuine long-range order implying that (bj) = 0. 
However, this form of mean-field approximation tacitly as- 
sumes that (bj) can be non-zero and so breaks the global 
U(l) symmetry of the BHM while preserving its transla- 
tional symmetry. In fact this complex scalar provides with 
a homogeneous order parameter <P = (bj) for all sites j 
accounting for the "mean-field" effect of hopping between 
neighbouring sites. This is formally identical to perform- 
ing a variational calculation using a translationally invari- 
ant Gutzwiller ansatz 02311221 for the ground state of the 
form 

i*i>=ni^>j' (4) 

3 




0.35 



Figure 2 (a) A block of L sites is represented together with 
its nearest neighbours from the two contiguous blocks. The 
hopping between the two sites at the end of the central 
block with the sites outside of it are approximated using 
a homogeneous SF order parameter ^ = (b\) = (b^) 
which is determined self-consistently 12411251 . In doing so, 
the Hamiltonian of the central block is decoupled from its 
neighbours sites rendering the system numerically solv- 
able, (b) Several MI lobes are shown which were computed 
using mean-field decoupling on the hopping between adja- 
cent ID blocks of L sites. The size of the block vary from 
the conventional L = 1 up to L = 8 site clusters. The 
solid gray lines are the boundaries from a 12th order SCP 
expansion which are included as a reference. 



over all sites j. Here | ipi) is the ground state of the single- 
site Hamiltonian Hi which can be determined straight- 
forwardly by numerical diagonalization once an appropri- 
ate occupancy cut-off is introduced [21 4|. By iterating 
this procedure a self-consistent solution is found where 
the <S> appearing in Hi is equal to the expectation value 
(ipi | bj | ipi) of its ground state and is equivalent to min- 
imizing the ground state energy ei(^) with respect to 
<P. The MI phase is identified by a symmetry preserving 
ground state | ipi) with ^ = and so its description is lim- 
ited to integer filled Fock states only, whilst the SF phase 
is described by a symmetry breaking superposition of on- 
site Fock states with <P ^ 0. Due to the product nature of 
the ansatz | ^i) no quantum correlations between sites are 
described by this approximation. 
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The rounded MI lobes predicted by the conventional 
single-site MF was shown earlier in Fig. 02a). We now ex- 
amine the results produced by a multi-site generalization of 
this approach I24II25I . Specifically, we apply the same de- 
coupling approximation to the hopping terms connecting 
the boundaries of each L-site blocks as shown schemati- 
cally in Fig. |2 a). A virtue of this approach is that the de- 
coupled systems now have a ID chain geometry and so the 
approximation departs from the effective infinite dimen- 
sional geometry of the single-site scheme. The MF decou- 
pling is equivalent to assuming a product ansatz 

i^)=ni^>c 

c 

over all L-site blocks C, with where | i/jl) is the ground 
state of the cluster Hamiltonian 

L 

Hl = -t($b\ + h.c.) - t(<Pb[ + h.c.) + 4i|<P| 2 - fi^hj 

3=1 

L-l L 
J=2 3=1 

Here the boundary MF is <P = (pi) = ipi) enforcing a 
homogeneous order parameter whose value is determined 
self-consistently as before. The MI and SF phase are again 
signalled by <P = and <P ^ 0, respectively. By parti- 
tioning the system in to clusters of L > 1 sites our de- 
scription loses its full translational symmetry, but in turn 
can describe MI and SF phases by more complex states 
with quantum correlations, such as particle-hole fluctua- 
tions, present within each block. 

As seen in Fig. |2jb) we find that the extension to a 
multi-site MF approach has only a marginal influence on 
both the shape and size of the predicted MI lobes. Going 
to large chains does extend the lobe to fractionally larger t, 
though the qualitative form of the MI lobes remain rounded 
similar to the single-site version. In principle as L — > oo 
this approach will converge to an exact diagonalization of 
the BHM. However, the results for numerically accessible 
chains sizes up to L = 8 show this convergence is ex- 
tremely slow and no re-entrant behaviour is visible. This 
suggests that the symmetry breaking induced by the bound- 
ary MF decoupling of the hopping is dominant and washes 
out vital physics underpinning re-entrance. 

3.2 Small system exact diagonalization Despite 
the fact that no finite-sized system possesses a genuine 
phase transition the exact diagonalization of small sys- 
tems can nonetheless provide important information about 
physics in the thermodynamic limit (26,27,28,29). Utiliz- 
ing symmetries and the sparseness of the BHM Hamilto- 
nian, with an on-site occupancy cut-off of 4 bosons and pe- 
riodic boundary conditions, we applied the Lanczos algo- 
rithm to numerically diagonalize system sizes from A = 2 
up to A = 11 sites [27). To examine the phase diagram we 
employed a quasi-grand-canonical approach. By working 




Figure 3 The sequence of particle and hole boundaries 
computed from small finite-sized clusters with periodic 
boundary conditions [ 26 , 29 ] . System sizes go from A = 2 
to A = 11 sites as the outermost to innermost lines, re- 
spectively. The solid gray lines are the boundaries for a 
12th order SCP expansion which are included as a refer- 
ence. 



in the canonical ensemble for an A site system we com- 
puted, as a function of t, the ground state energy ejv+i(i), 
ejv(i) and ejv-i(f) for the system with exactly A + 1, A 
and A — 1 bosons, respectively. We then extract an effec- 
tive grand-canonical ensemble phase diagram by comput- 
ing the particle and hole boundaries in the (t, /i) plane as 
the degeneracy curves between these particle number sec- 
tors given by the solutions to the equations 

ejv+i(i) - e N (t) + [i = 0, (particle), 
CAr(*)-ejv-i(*)+A» = 0. (hole) 

In Fig.[3]the resulting particle and hole boundaries for the 
MI lobe for increasing system sizes A are shown. An in- 
dication that these curves do not produce genuine phase 
boundaries is that the particle and hole boundaries do not 
intersect for any finite A and so predict a MI lobe of in- 
finite size. However, a key feature, clearly visible for all 
system sizes examined, is that the hole boundary already 
shows concavity at moderate values of t characteristic of 
re-entrance. Moreover there is a visible tendency of the 
finite-size system boundaries to slowly converge to the ex- 
pected MI lobe in the thermodynamic limit as predicted by 
the 12th order SCP expansion. This feature is the basis of 
finite-size scaling calculations |26 j29| . 

Both the finite-sized cluster and the multi-site MF ap- 
proximations are based on an ansatz with a repeated clus- 
ter of variable size which breaks the underlying transla- 
tional symmetry of the BHM Hamiltonian in Eq. (Q]). Yet 
the phase boundaries predicted for each of these meth- 
ods are rather different. We see that the finite-sized cluster 
calculation produces systematic overestimations of the MI 
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lobe for any cluster size while the multi-site MF approxi- 
mation produces a systematic underestimation. The reason 
for this is that unlike MF decoupling the finite-sized clus- 
ter approach does not break the U(l) global symmetry of 
the BHM and is unable to mimic the quasi-long range or- 
der characteristic of a SF thereby favouring the MI phase. 
In contrast the multi-site MF approximation can break this 
symmetry favouring the SF phase. That a precursor to re- 
entrance is visible in finite-sized clusters demonstrates that 
a symmetry preserving description of kinetic energy is cru- 
cial to its emergence in an approximation. This conclusion 
is further enhanced by results obtained with more sophis- 
ticated approaches such as the variational cluster pertur- 
bation method lf30l . There the system is partitioned into 
finite-size clusters, like that used in the multi-site MF ap- 
proximation, however hopping between the clusters is in- 
stead treated perturbatively. Crucially this allows the U(l) 
symmetry to be preserved in both the SF and MI phases. 
This phase diagram produced by this method then resem- 
ble the exact-diagonalization results here and also similarly 
predict re-entrance for any cluster size ll30ll . 

3.3 Real-space renormalization group The expo- 
nential growth in the Hilbert space with N means that exact 
diagonalization cannot be feasibly pushed much beyond 
the sizes considered here already. For this reason we ex- 
amine the predictions of a real-space RG approach for the 
BHM following the scheme originally introduced by Singh 
and Rokhsar OTl some time ago. In that work they focused 
on the canonical ensemble at unit-filling to analyse the KT 
transition point at the tip. Here for the first time we refor- 
mulate their scheme within the quasi-grand-canonical en- 
semble approach, just applied to the finite-sized clusters, 
to obtain an RG prediction of the complete phase diagram 
of the BHM. The RG approximation combines two crucial 
features of the MF and finite-sized clusters considered al- 
ready. First, like finite-sized clusters it preserves the U(l) 
global symmetry of the BHM. Second, like the MF decou- 
pling RG scheme curtails the exponential growth of the 
system's Hilbert space enabling the thermodynamic limit 
to be reached. 

The real-space RG scheme applied is based on par- 
titioning the ID system into non-overlapping two site 
blocks. In contrast to the other approximations the RG 
scheme OTl simplifies the BHM by truncating the on- 
site Fock basis to the occupation states | 0) , 1 1) and | 2) 
thereby excluding triple or higher occupancies. In this ba- 
sis the kinetic contribution K to the block Hamiltonian 
Hb, describing hopping between adjacent sites within a 
block, is 

£ = -^(|01}(10|+h.c.) 

-i o (|ll)(20| + |ll) (02|+h.c.) 

-ip(|21><12|+h.c.) 

where th and t p , specify the hopping amplitude for a hole 
or additional particle, respectively, while t a is the ampli- 
tude for the creation or annihilation of a particle-hole pair. 



When constructing K directly from the hopping terms in 
the BHM in Eq. (Q} we have that t h = t /\f2 = t p /2 = t 
reflecting the fact that the BHM is not particle-hole sym- 
metric. The block Hamiltonian Hi, is then completed by 
adding the interaction and chemical potential terms, which 
in the truncated basis are diagonal on-site operators with 
diagonal matrix elements (0, 0, U) and (0, /i, 2/x), respec- 
tively EQ. 

The RG scheme begins by diagonalizing Hi, and identi- 
fying three eigenstates | 0') , 1 1') and | 2') which represent 
the lowest energy state with one, two and three bosons per 
block, respectively. The Fock state truncation enables these 
states to be readily determined as 

|0'> = -^(|01> + |10», 

|l') = cos(fl)|ll> + -^sin(fl)(|02) + |20», 

|2') = ^(U2> + |21», 

with corresponding energies 

Eo> = —fJ> — th, 

E v = -2(jl- \/2 to sin(20) + Usm 2 (9), 
E 2 > = -3/x -tp + U, 

where 6 = ^tsai" 1 (2^/2t /U). Using this appropriately 
chosen reduced basis of states for a single block we can 
now interpret the block as a new super-site and formulate 
a renormalized block Hamiltonian H' b spanning two such 
adjacent super-sites. In this way the RG scheme dramati- 
cally limits the degrees of freedom retained. The Hamilto- 
nian H' b has exactly the same type of terms as Hi,, once the 
replacements | 0) — > | 0') , etc. have been made, but pos- 
sesses different matrix elements t' h ,t' ,t' , U' and //. The 
new on-site matrix elements are given by 

// = E v -E , = fi- t h [l - 2sin(20)] + t/sin 2 (6>), 
U' = E 2 > + E& - 2E V = Uscc(29) -t p - t h . 

The new hopping amplitudes are then determined by com- 
puting directly in the two block Hilbert space the matrix 
elements 

t h = -(Q?l'\k\l'Q?), 

t' a = -( 1' 1' 2' 0') = -(1'1'|/?|0'2'), 

t' p = -(2'l'\K\l'2'), 

where K is the hopping term connecting the two blocks RTl . 
Repeating this process for k RG iterations yields ma- 
trix elements ,*o ,t p ,U^ k \ fj,^ as well as en- 
ergies Eq , Ei , E^ and their corresponding states 
|0W),|lW),|2«). 

Typically after k w 100 RG iterations the scheme 

(k) 

drives the system to particle-hole symmetry with t h = 
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(a) t' t' (c) 




Figure 4 (a) The real-space RG scheme adopted here J3TJ. The Hamiltonian of a two site block is diagonalized and the 
lowest one, two and three boson eigenstates are kept. This forms an isometry mapping 3 <g> 3 — > 3 which renormalizes 
a block of two sites into a single super-site coupled to it neighbours with a different hopping, i.e. t — > t'. (b) This 
RG scheme is equivalent to variationally minimizing a TTN ansatz. Here each isometry is represented as a 3 leg tensor 
(triangle) composed as a hierachical network, (c) The MI lobe predicted by the real-space RG scheme is displayed, along 
with the green point marking the predicted KT critical point [31]. The inset is a zoom of phase boundaries around the 
tip highlighting the observed re-entrance. The solid gray lines are again the boundaries from a 12th order SCP expansion 
included as a reference. 



to = tp = t 1 -^ and converges to one of two fixed 
points. The MI is signified by (t/U)^ = whereupon 
| 1«) = | l( fc -!) 1 (*-!)) illustrating that the RG has re- 
moved the short-range correlations and converged to a 
unit-filled Fock state. The SF is signified by the ratio 
(t/J7)'*' = 0.5165 possessing non-zero particle-hole fluc- 
tuations. The repeated renormalization of two-site blocks 
is equivalent to variationally approximating the ground 
state of the BHM via a hierarchically structured tensor tree 
network (TTN), as schematically shown in Fig. H|b). As 
with the multi-site MF decoupling and finite-sized clusters 
this TTN ansatz breaks the translationally symmetry of the 
BHM. 

The RG ansatz possesses crucial similarities to a finite- 
sized cluster approach. If after k RG steps the single open 
leg at the top of the TTN was capped by the final block state 
1 1W) the resulting network is a highly efficient represen- 
tation of a many-body state of an N = 2 k site chain filled 
with exactly N bosons. It therefore preserves the U(l) 
global symmetry. Similarly capping the TTN with | 0*- fe ^) 
or | 2^) yields a representation of a state with N — 1 and 
N + 1 bosons, respectively. Thus, while the RG procedure 
can easily reach the thermodynamic limit it can only de- 
scribe very limited ±1 particle-number fluctuations due to 
the repeated truncation applied. This means that the RG ap- 
proach might capture correlations of a MI reasonably, but 
can only offer a very poor description of the SF phase. 

In an analogous way to the finite-size cluster approxi- 
mation, the RG prediction for the phase diagram was deter- 
mining by numerically solving for the degeneracy curves 



E W_ E W = Q and E W_ E W _ gin the (t, fi) plane be- 
tween | l( fe )) and the states | Q( k >) and | 2W), respectively. 
The resulting hole and particle boundaries are reported in 
Fig. E|c). The point where they meet forms the lobe tip 
whose location is a unstable fixed point that can be found 
analytically to be at (t/U)^ = 0.215 (3TJ- The predicted 
MI lobe has a triangular shape and from inset of Fig. Etc), 
zoomed in near the tip, the appearance of re-entrance is 
evident. Compared to the SCP expansion result the parti- 
cle boundary is in reasonable agreement, however the KT 
point is under-estimated and the hole boundary differs no- 
ticeably away from the f 4 limit. Nonetheless the RG 
approximation has captured the re-entrant feature of the 
MI lobes in ID and its success strongly points to a ?7(1) 
symmetric description of particle-hole fluctuations of a MI 
state as being essential. This conclusion is consistent with 
recent investigations linking re-entrance in the BHM in ID 
to entanglement using a more sophisticated infinite matrix 
product ansatz ll35ll . 

4 Conclusions Re-entrance is often difficult to cap- 
ture in approximations of the ID BHM. Nevertheless, we 
have shown that this real-space RG scheme does in fact re- 
produce this phenomenon despite being much simpler than 
the DMRG, high order SCP expansions or QMC. The fail- 
ures of a multi-site MF decoupling and finite-sized cluster 
approach reveal clues as to why the RG approximation has 
this ability. Specifically, it retains some, albeit limited, ca- 
pacity to describe in a particle-number symmetric way rel- 
evant fluctuations of the MI while truncating sufficiently 
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to enable the thermodynamic limit to be approached. This, 
and the precursors to re-entrance seen in finite-sized clus- 
ters, highlights the importance of short range correlations 
in the MI state near the tip. Furthermore, the truncation 
of the on-site occupancy in the RG scheme shows that 
particle-hole pair correlations in particular, arising from 
kinetic energy, play a significant role. This can be con- 
trasted with the related issue of the location of the KT criti- 
cal point where long-range correlations are instead crucial. 
An interesting open question is how an extended real-space 
RG scheme, where the initial on-site Fock basis and subse- 
quent RG steps truncate less, improves upon the simplest 
approach used here while retaining a physically intuitive 
form. Our work here suggests that a moderate improve- 
ment in the description of particle-hole fluctuations could 
substantially improve the predicted MI lobe shape and KT 
point. 
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